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Spectral Functions from Hadronic r Decays and QCD 
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Hadronic decays of the r lepton provide a clean environment to study hadron dynamics in an energy regime 
dominated by resonances, with the interesting information captured in the spectral functions. Recent results from 
ALEPH on exclusive channels are presented, with some emphasis on the nir° final state which plays a crucial role 
for the determination of the hadronic contribution to the muon anomalous magnetic moment. Inclusive spectral 
functions are the basis for QCD analyses, delivering an accurate determination of the strong coupling constant, 
quantitative information on nonperturbative contributions and a measurement of the mass of the strange quark. 



1. Introduction 

Hadrons produced in r decays are born out of 
the charged weak current, i.e. out of the QCD 
vacuum. This property guarantees that hadronic 
physics factorizes in these processes which are 
then completely characterized for each decay 
channel by spectral functions as far as the total 
rate is concerned . Furthermore, the produced 
hadronic systems have 1 = 1 and spin-parity 
J p = 0+, I" (V) and J p = 0~, 1+ (A). The spec- 
tral functions are directly related to the invariant 
mass spectra of the hadronic final states, normal- 
ized to their respective branching ratios and cor- 
rected for the r decay kinematics. For a given 
spin-1 vector decay, one has 
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where V ud = 0.9748 ± 0.0010 denotes the CKM 
weak mixing matrix element [1] and Sew ac- 
counts for electroweak radiative corrections (see 
below). Isospin symmetry (CVC) connects the 
r and e + e~ annihilation spectral functions, the 
latter being proportional to the R ratio. 

Hadronic r decays are a clean probe of hadron 
dynamics in an interesting energy region dom- 



inated by resonances. However, perturbative 
QCD can be seriously considered due to the rel- 
atively large r mass. Many hadronic modes have 
been measured and studied, while some earlier 
discrepancies (before 1990) have been resolved 
with high-statistics and low-systematics experi- 
ments. Samples of ~ 4 x 10 5 measured decays 
are available in each LEP experiment and CLEO. 
Conditions for low systematic uncertainties are 
particularly well met at LEP: measured samples 
have small non-r backgrounds (~ 1%) and large 
selection efficiency (92%), for example in ALEPH. 

Recent results in the field are discussed in this 
report. 

2. New ALEPH Spectral Functions 

Preliminary spectral functions based on the full 
LEP1 statistics are available from ALEPH. The 
corresponding results for the branching fractions 
are given separately [2]. The analysis uses an 
improved treatment of photons as compared to 
the published analyses based on a reduced sam- 
ple [3,4]. Spectral functions are unfolded from 
the measured mass spectra after background sub- 
traction using a mass-migration matrix obtained 
from the simulation in order to account for detec- 
tor and reconstruction biases. Backgrounds from 
non-rr events are small (< 1% and subtractions 
are dominated by r decay feedthroughs. This is 



2 



- MC+z-bg (all, 
z-bg (all) 
-,(e) 



I 10 ' 



3 ^ 


T-/>£ (71) 








/ - 


1 , i , , 1 , 


, , 1 , , 




Mass 2 (Gey 2 ) 



Figure 1. Feedthrough in the tttt channel from 
other r decay modes in the ALEPH analysis, 
as estimated from the Monte Carlo simulation. 
The light-shaded histogram indicates the total 
feedthrough, while the dark-shaded histograms in 
each plot show the dominant individual sources. 



illustrated in Fig. 1 for the tttt channel, where 
the main feedthrough contributions are from the 
7r mode at large masses (extra photon from radia- 
tion and fake photons from hadronic interactions 
in the electromagnetic calorimeter) and mostly 
the a\ — > 7r27r° mode (one tt° lost). 

The raw and unfolded spectra for the tt~tt° 
mode are given in Fig. 2. Apart from the normal- 
ization effects studied in the branching ratio anal- 
ysis [2] , systematic biases affecting the mass spec- 
tra, are studied separately. The corresponding 
uncertainties and correlations in mass are incor- 
porated in covariance matrices for every consid- 
ered channel. According to the number of pions 
in the final state, the spectral functions are sep- 
arated into vector and axial-vector components. 
As for final states involving KK pairs, specific in- 
put is required to achieve the V~ A separation [5]. 
Strange final states are treated separately. 
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Figure 2. Spectral functions in the tttt 
mode from ALEPH before (histogram) and after 
(points) the unfolding procedure. 



3. Specific Final States 

3.1. The 2?r Vector State 

3.1.1. The Data 

The spectral function from r — > is t tt~tt° in the 
full-LEPl ALEPH analysis (~ 10 5 events) is pre- 
sented in Fig 3. It is dominated by the p(770) 
resonance, with some broad structure around 
1.25 GeV. 

Results from CLEO [6] are in good agreement 
with the ALEPH data. The statistics is compara- 
ble in the two cases, however whereas the accep- 
tance is very flat across the mass spectrum in the 
case of ALEPH, a strong increase is observed for 
the CLEO case. As a consequence, ALEPH data 
are more precise below the p peak, while CLEO is 
more precise above. Thus the two pieces of data 
are complementary. The comparison is given in 
Fig. 4. 

3.1.2. tttt Spectral Functions and tt Form 
Factors 

It is useful to carefully write down all the fac- 
tors involved in the comparison of e + e~ and r 
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Figure 3. The pion form factor squared |F~(s)| 2 
measured by ALEPH. Because of the unfolding 
procedure, the errors in nearby bins are strongly 
correlated (only diagonal errors are shown). 



spectral functions in order to make explicit the 
possible sources of CVC breaking. On the e + e~ 
side we have 
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where /3g(s) is the threshold kinematic factor and 
F®(s) the pion form factor. On the r side, the 
physics is contained in the hadronic mass dis- 
tribution and one can write correspondingly to 
Eq. (3) 
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Figure 4. The comparison between ALEPH and 
CLEO spectral functions in the irir mode. 



C(s) = 

SU(2) symmetry (CVC) implies 
v-(s) = v (s) 



2s 
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It is important to note that experiments on 
r decays measure the rate inclusive of radiative 
photons, i.e. for r — > v r irTT° (j) . The measured 
spectral function is thus v*_(s) = v-(s) G(s), 
where G(s) is a radiative correction. Also, as 
|F"(s)| 2 should only involve the pion hadronic 
structure, it must not contain vacuum polariza- 
tion in the photon propagator. 

Three levels of SU(2) breaking can be identi- 
fied: 

• electroweak radiative corrections to t decays 
are contained in the Sew factor [7,8] which 
is dominated by short-distance effects. As 
such it is expected to be weakly depen- 
dent on the specific hadronic final state, 
as verified by detailed computations in the 
t — ► (tt,K)i/ t channels [9]. Recently, de- 
tailed calculations have been performed for 
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the 7T7r° channel [10] which also confirm the 
relative smallness of the long-distance con- 
tributions. One can write the total correc- 
tion as 



chad chad 
C _ °EW°EM 
b EW ~ -JiJ 

'-'em 



(5) 



where S E ^ is the leading-log short-distance 
electroweak factor (which vanishes for lep- 
tons) and S } ^ I ' lcp are the nonleading elec- 
tromagnetic corrections. The latter correc- 
tions are calculated in Ref. [8] at the quark 
level and in Ref. [10] at the hadron level for 
the 7T7r° decay mode, and in Ref. [7,8] for 
leptons. The total correction amounts [13] 
to S^ u = 1.0198 ±0.0006 for the inclusive 
hadron decay rate and = (1.0232 ± 

0.0006) S%*m(s) for the tttt decay mode, 
where S E \ { (s) is an s-dependent radiative 
correction [10]. 

• pion mass splitting, which is almost com- 
pletely from electromagnetic origin, directly 
breaks isospin symmetry in the spectral 
functions [11,12] since /3_(s) ^ (3o(s). 

• symmetry is also broken in the pion form 
factor [11,10]. The p width is affected by 7r 
and p mass splittings and by explicit elec- 
tromagnetic decays such as irj, 777, 

and 7T7T7. Isospin violation in the strong 
amplitude is expected to be negligible be- 
cause of the small absolute mass difference 
between u and d quarks. 

3.1.3. Fitting the 27r Spectral Function 

The 2ir spectral function is dominated by the 
p resonance. This provides us with an opportu- 
nity to study the relevant description of a wide 
hadronic resonance. The Gounaris-Sakurai [14] 
(GS) parametrization satisfies analyticity and 
unitarity through finite width corrections. The 
line shape has to be modified to account for the 
effect of p — uj interference. 

The pion form factor is fitted with interfering 
amplitudes from p, p' and p" vector mesons with 
relative strengths 1, (3 and 7. A phase <j)p is also 



considered, since the relative phase of the p and 
p' amplitude is a priori unknown. In practice we 
fit the F®(s) form factor from e + e~ data and the 
F~(s) form factor from the r spectral function 
duly corrected for SU(2) breaking, however only 
in the spectral function phase space and in the 
Sew factor. In this way, the two form factors can 
be readily compared, i.e. the masses and widths 
of the dominating p resonance in the two isospin 
states can be unambiguously determined. On the 
e + e~ side, all available data have been used (com- 
plete references can be found in Ref. [13], except 
those of NA7 [15] which are questionable [16]. 
However, in the p mass region the new precise 
results from CMD-2 [17] are dominating. On the 
r side, the accurate data from ALEPH and CLEO 
are used. 

The systematic uncertainties are included in 
the fits through appropriate covariance matrices. 
The p mass systematic uncertainty in t data is 
mostly from calibration (0.7 McV for ALEPH 
and 0.9 MeV for CLEO). The corresponding un- 
certainties on the p width are 0.8 and 0.7 MeV. 
Both mass and width determinations are system- 
atically limited in r data, however they are uncor- 
rected between ALEPH and CLEO. Due to the 
large event statistics, the fits are quite sensitive to 
the precise line shape and to the interference be- 
tween the different amplitudes. Of course, p — u 
interference is included for e + e~ data only and 
the corresponding amplitude (a puJ ) fitted. Fits 
are performed separately for the r and e + e~ form 
factors. The upper range of the fit is taken at 3.6 
GeV 2 for e+e", but at 2.4 GcV 2 for r to avoid 
the extreme edge of the spectral function where 
large corrections are applied. 

It turns out that the fitted p masses and widths 
(essentially the latter) are quite sensitive to the 
strength of the p' and p" amplitudes, (3 and 7. So, 
depending on the type of fit, the derived values 
exhibit some systematic shifts. The e+e~ and 
t fits yield significantly different p' amplitudes 
and phases: the phase in the r fit is consistent 
with 180°, while it comes out 30° smaller in the 
e + e~ fit. This is merely a reflection of the dis- 
crepancy between the measurements of the pion 
form factor, especially between 0.8 and 1.0 GeV. 
A fit to both data sets is also performed, keep- 
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Figure 5. The pion form factor squared |F°(s)| 2 Figure 6. The pion form factor squared \F n (s)\ 2 

from e + e~ data and the combined fit from Ta- from r data and the combined fit from Table 1. 

ble 1. 



ing common values for the p' and p" parameters 
(second-order with respect to the dominant in- 
vestigated p parameters). Table 1 presents the 
results of the common fits, the quality of which 
can be inspected in Figs. 5 and 6. 

The differences found between the masses and 
widths of the charged and neutral p's are some- 
what larger than predicted by Chiral Perturba- 
tion Theory (%PT) [18](|m p - - m p o \ < 0.7 MeV) 
and isospin breaking (T p - — r p o = (0.7 ± 0.3) 
MeV). However, if the mass difference (2.6 MeV) 
is taken as an experimental fact, then the larger 
observed width difference (3.1 MeV) could be ex- 
plained by xPT (0.7 MeV from isospin breaking 
+ 1.5 MeV from the mass difference). Whether 
or not these discrepancies are of experimental or 
theoretical nature is still open [13,19] and needs 
to be further investigated. 



3.2. The 4tt Vector State 

The 47T final states have also been studied [3, 
20]. Tests of CVC are severely hampered by 
large deviations between different e + e~ experi- 
ments which disagree beyond their quoted sys- 
tematic uncertainties. A CLEO analysis studies 
the resonant structure in the 3irir° channel which 
is shown to be dominated by unr and a\ir contri- 
butions. The lutt spectral function is interpreted 
by a sum of p-like amplitudes. The mass of the 
second state is however found at (1523 ± 10)MeV, 
in contrast with the value (1406 ± 14)MeV from 
the fit of the 2tt spectral function. This point 
has to be clarified. Following a limit of 8.6% ob- 
tained earlier by ALEPH [21], CLEO sets a new 
95% CL limit of 6.4% for the relative contribution 
of second-class currents in the decay r — > v T ix~ u 
from the hadronic angular decay distribution. 

The new ALEPH spectral function in the 3irir° 
mode is shown in Fig. 7: it compares well with 
the CLEO published results [20], except maybe 
in the threshold region near 1 GeV. 
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773.3 ± 0.6 


772.7 ± 0.5 


T p o 




145.2 ± 1.3 


146.4 ± 0.9 


m p - 


775.0 ± 0.6 




775.3 ± 0.6 




149.5 ± 1.1 




149.5 ± 0.8 






(2.02 ±0.10) 10~ 3 


(1.98 ±0.10) 10- 3 


a 


0.195 ± 0.028 


0.123 ± 0.011 


0.172 ± 0.006 




173.0 ± 7.0 


139.4 ± 6.5 


178.2 ± 4.5 


m p > 


1440 ± 34 


1337 ± 35 


1415 ± 15 


I> 


597 ± 102 


569 ± 81 


528 ± 42 
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0.095 ± 0.029 


0.048 ± 0.008 


0.072 ± 0.006 




0. 


0. 


0. 


m p >> 


1713 


1713 ± 15 


1741 ± 20 




235 


235 


235 



Table 1 

Results of fits to the pion form factor squared to r and e + e~ data (ALEPH and CLEO) separately, then 
combined. The parametrization of the p line shape (Gounaris-Sakurai) is described in the text. All mass 
and width values are in MeV and the phase (ftp in degrees. 



The r spectral function in the four-pion modes 
can be used to predict the corresponding isospin- 
rotated e + e~ cross sections using the SU(2) sym- 
metric relations 

1=1 o 4?ra2 

/=! 47TO 2 
"e+e--* tt+tt-tt°it ~ — ~ l v 2n - IT + TT v it- 3tt° 

Comparisons can be performed after applying 
known isospin-breaking corrections [12,13]. As 
shown in Figs. 8 and 9, agreement is observed in 
the 2ir + 2ir~ mode, whereas the situation is con- 
fused in the 7r + 7r~27r° mode where a large spread 
exists between the different e + e~ results (refer- 
ences can be found in Ref. [13]). 

3.3. Axial-vector States 

The decay r — ► v t 3tt is the cleanest place to 
study axial- vector resonance structure. The spec- 
trum is dominated by the 1 + a\ state, known to 
decay essentially through p-n. A comprehensive 
analysis of the tt~2tt° channel has been presented 
by CLEO. First, a model- independent determi- 
nation of the hadronic structure functions gave 
no evidence for non-axial- vector contributions (< 
17% at 90% CL) [22]. Second, a partial- wave am- 



plitude analysis was performed [23]: while the 
dominant pir mode was of course confirmed, it 
came as a surprize that an important contribu- 
tion (~ 20%) from scalars (a, /o(1470), / 2 (1270)) 
was found in the 2-7T system. 

The ai — ► 7r~27T° precisely determined line 
shape shows the opening of the K*K decay mode 
in the total a\ width. The derived branching ra- 
tio, B( ai -» K*K) = (3.3 ± 0.5)% is in good 
agreement with ALEPH results on the KKir 
modes which were indeed shown (with the help 
of e + e~ data and CVC) to be axial-vector (ai) 
dominated with B(ai -» K*K) = (2.6 ± 0.3)% 
[5] . No conclusive evidence for a higher mass state 
(a'j) was found in this analysis. 

4. Inclusive Nonstrange Spectral Func- 
tions 

The r nonstrange spectral functions have been 
measured by ALEPH [3,4] and OPAL [24]. The 
procedure requires a careful separation of vector 
(V) and axial-vector (A) states involving the re- 
construction of multi-7r° decays and the proper 
treatment of final states with a K K pair. The 
V and A spectral functions are given in Fig. 10. 
They show a strong resonant behaviour, domi- 
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Figure 7. The spectral function for the Sirn' 
mode as determined by ALEPH and CLEO. 



Figure 8. The cross section for e + e _ — > 2ir + 2ir~ 
from e + e~ experiments and derived from the 
7r37T° spectral function measured by ALEPH. 



5. QCD Analysis of Nonstrange r Decays 



nated by the lowest p and a\ states, with a ten- 
dancy to converge at large mass toward a value 
near the parton model expectation. Yet, the vec- 
tor part stays clearly above while the axial-vector 
one lies below. Thus, the two spectral functions 
are clearly not 'asymptotic' at the r mass scale. 

The new inclusive ALEPH V and A spectral 
functions are given in Figs. 11 and 12 with a 
breakdown of the respective contributions. 

The V + A spectral function, shown in Fig. 13 
has a clear pattern converging toward a value 
above the parton level as expected in QCD. In 
fact, it displays a textbook example of global du- 
ality, since the resonance-dominated low-mass re- 
gion shows an oscillatory behaviour around the 
asymptotic QCD expectation, assumed to be 
valid in a local sense only for large masses. This 
feature will be quantitatively discussed in the 
next section. 



5.1. Motivation 

The total hadronic r width, properly normal- 
ized to the known leptonic width, 



r(r- 



hadrons v T ) 



r( 7 



V e V T 



(6) 



should be well predicted by QCD as it is an inclu- 
sive observable. Compared to the similar quan- 
tity defined in e + e~ annihilation, it is even twice 
inclusive: not only are all produced hadronic 
states at a given mass summed over, but an inte- 
gration is performed over all the possible masses 
from mjr to m T . 

This favourable situation could be spoiled by 
the fact that the Q 2 scale is rather small, so that 
questions about the validity of a perturbative ap- 
proach can be raised. At least two levels are to be 
considered: the convergence of the perturbative 
expansion and the control of the nonperturbative 
contributions. Happy circumstances make these 
latter components indeed very small [25,26]. 
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Figure 9. The cross section for e + e _ — > ir + ir~2iT 
from e + e~ experiments and derived from the 
7r37r° and 37T7r° spectral functions measured by 
ALEPH. 



5.2. Theoretical Prediction for R T 

The imaginary parts of the vector and 
axial-vector two-point correlation functions 

n4d' V/a( s )' w ^ ^ e s P^ n ^ °^ ^he nacu " on i c 
system, are proportional to the r hadronic spec- 
tral functions with corresponding quantum num- 
bers. The non-strange ratio R T can be written 
as an integral of these spectral functions over 
the invariant mass-squared s of the final state 
hadrons [27]: 

so 2 

R T (s ) = 12^5 EW f-(l--) (7) 
J s \ s J 



1 + 2—1 Imn (1) (s + ie) + ImII (0) (s + ie) 
so 



By Cauchy's theorem the imaginary part of 

is proportional to the discontinuity across the 

positive real axis. 

The energy scale s = is large enough so 
that contributions from nonperturbative effects 
are small. It is therefore assumed that one can 
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Figure 10. Inclusive nonstrange vector (top) 
and axial- vector (bottom) spectral functions from 
ALEPH and OPAL. The dashed line is the expec- 
tation from the naive parton model. 
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Figure 11. Preliminary ALEPH inclusive vector 
spectral function with its different contributions. 
The dashed line is the expectation from the naive 
parton model. 



Figure 12. Preliminary ALEPH inclusive axial- 
vector spectral function with its different contri- 
butions. The dashed line is the expectation from 
the naive parton model. 



use the Operator Product Expansion (OPE) to 
organize perturbative and nonperturbative con- 
tributions [28] to R t (sq)- 

The theoretical prediction of the vector and 
axial-vector ratio R T .v/A can thus be written as: 
3 

Rt.V/A = -^\Vud\ ^ EW (®) 

x ( i + * (0) + aw + €7/7" + E Cv a) 

\ D=4,6,8 J 

with the residual non-logarithmic electroweak 
correction S' EW = 0.0010 [8], neglected in the fol- 
lowing, and the dimension D = 2 contribution 
^ud vJa^ f rom quark masses which is lower than 
0.1% for u, d quarks. The term 5^°' is the purely 
perturbative contribution, while the 5^ are the 

— D/2 

OPE terms in powers of s of the following 
form 

AD) {Oud)y/A , s 

°ud,V/A (_ )D/2 

dimO=D y V > 

where the long-distance nonperturbative effects 



are absorbed into the vacuum expectation ele- 
ments (O u d). 

The perturbative expansion (FOPT) is known 
to third order [29]. A resummation of all known 
higher order logarithmic integrals improves the 
convergence of the perturbative series (contour- 
improved method FOPTci) [30]. As some ambi- 
guity persists, the results are given as an average 
of the two methods with the difference taken as 
a systematic uncertainty. 

5.3. Measurements 

The QCD analysis of the r hadronic width has 
not yet been completed with the final ALEPH 
spectral functions. Results given below corre- 
spond to the published analyses with a smaller 
data set. 

The ratio R T is obtained from measurements of 
the leptonic branching ratios: 

R T = 3. 647 ±0.014 (10) 

using a value B{t~ — > v e v T ) — (17.794 ± 
0.045)% which includes the improvement in ac- 
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Figure 13. The preliminary inclusive V + A 
nonstrange spectral function from ALEPH. The 
dashed line is the expectation from the naive par- 
ton model. 



curacy provided by the universality assumption 
of leptonic currents together with the measure- 
ments of J5(r _ — > e~ v e v T ), B(t~ — > fi~ v^v^) 
and the r lifetime. The nonstrange part of R T is 
obtained by subtracting out the measured strange 
contribution (see last section). 

Two complete analyses of the V and A 
parts have been performed by ALEPH [4] and 
OPAL [24]. Both use the world-average leptonic 
branching ratios, but their own measured spec- 
tral functions. The results on a s (m 2 .) are there- 
fore strongly correlated and indeed agree when 
the same theoretical prescriptions are used. 

5.4. Results of the Fits 

The results of the fits are given in Table 2 for 
the ALEPH analysis. Similar results are obtained 
by OPAL. It is worth emphasizing that the non- 
perturbative contributions are found to be very 
small, as expected. The limited number of ob- 
servables and the strong correlations between the 
spectral moments introduce large correlations, es- 
pecially between the fitted nonperturbative oper- 



ALEPH 


a. s (m 2 T ) 


5np 


V 
A 


0.330 ±0.014 ±0.018 
0.339 ±0.013 ±0.018 


0.020 ± 0.004 
-0.027 ± 0.004 


V+A 


0.334 ± 0.007 ± 0.021 


-0.003 ± 0.004 



Table 2 

Fit results of a s (m^) and the OPE nonperturba- 
tive contributions from vector, axial-vector and 
(V + A) combined fits using the corresponding 
ratios R T and the spectral moments as input pa- 
rameters. The second error is given for theoretical 
uncertainty. 



ators. 

One notices a remarkable agreement within sta- 
tistical errors between the a s (m^.) values using 
vector and axial-vector data. The total nonper- 
turbative power contribution to R T ^y+A is com- 
patible with zero within an uncertainty of 0.4%, 
that is much smaller than the error arising from 
the perturbative term. This cancellation of the 
nonperturbative terms increases the confidence 
on the a s (m^) determination from the inclusive 
(V + A) observables. 

The final result from ALEPH is : 

a s {m 2 T ) = 0.334 ±0.007 oxp ± 0.021th (11) 

where the first error accounts for the experimen- 
tal uncertainty and the second gives the uncer- 
tainty of the theoretical prediction of R T and the 
spectral moments as well as the ambiguity of the 
theoretical approaches employed. 

In the OPAL analysis the corresponding re- 
sults arc quoted within three prescriptions for the 
perturbative expansion, respectively FOPTci, 
FOPT and with renormalon chains, i.e. 

a s {m 2 T ) = 0.348 ±0.009 oxp ± 0.019th (12) 
= 0.324 ±0.006 oxp ± 0.013th (13) 
= 0.306 ±0.005 oxp ± 0.011th (14) 

5.5. Test of the Running of a s (s) at Low 
Energies 

Using the spectral functions, one can simulate 
the physics of a hypothetical r lepton with a 
mass y/so smaller than m T through equation (7) 
and hence further investigate QCD phenomena 
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Figure 14. The ratio i? T ^ +J 4 versus the square "r 
mass" so. The curves are plotted as error bands 
to emphasize their strong point-to-point correla- 
tions in so. Also shown is the theoretical predic- 
tion using the results of the fit at s = to 2 . 



at low energies. Assuming quark-hadron duality, 
the evolution of R t (sq) provides a direct test of 
the running of a s (so), governed by the RGE (3- 
function. On the other hand, it is a test of the 
validity of the OPE approach in r decays. 
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Figure 15. The running of a s (so) obtained from 
the fit of the theoretical prediction to R Ti v + a{so)- 
The shaded band shows the data including ex- 
perimental errors. The curves give the four-loop 
RGE evolution for two and three flavours. 



5.6. Discussion on the Determination of 

a s {ml) 

The evolution of the a s (m 2 ) measurement from 
the inclusive (V + A) observables based on the 
Runge-Kutta integration of the differential equa- 
tion of the renormalization group to N 3 LO [31,33] 
yields for the ALEPH analysis 

a s (M|) = 0.1202±0.0008 e xp±0.0024th±0.0010ovoi(15) 



The functional dependence of R t ,v + a{so) is 
plotted in Fig. 14 together with the theoretical 
prediction using the results of Table 2. Below 
1 GeV 2 the error of the theoretical prediction of 
Rt,v+a{so) starts to blow up due to the increas- 
ing uncertainty from the unknown fourth-order 
perturbative term. Fig. 15 has the same phys- 
ical content as Fig. 14, but translated into the 
running of a s (so), i.e., the experimental value for 
a s (so) has been individually determined at ev- 
ery so from the comparison of data and theory. 
Good agreement is observed with the four-loop 
RGE evolution using three quark flavours. 

The experimental fact that the nonperturba- 
tive contributions cancel over the whole range 
1.2 GeV 2 < s < m 2 leads to confidence that 
the a s determination from the inclusive (V + A) 
data is robust. 



where the last error stands for possible ambigu- 
ities in the evolution due to uncertainties in the 
matching scales of the quark thresholds [33] . 

The result (15) can be compared to the pre- 
cise determination from the measurement of the 
Z width, as obtained in the global electroweak 
fit. The variable Rz has similar advantages to 
R T , but it differs concerning the convergence of 
the perturbative expansion because of the much 
larger scale. It turns out that this determination 
is dominated by experimental errors with very 
small theoretical uncertainties, i.e. the reverse of 
the situation encountered in r decays. The most 
recent value [34] yields a s (M|) = 0.1183±0.0027, 
in excellent agreement with (15). Fig. 16 illus- 
trates well the agreement between the evolution 
of a s (m 2 T ) predicted by QCD and a s (M^). 
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Figure 16. Evolution of the strong coupling (mea- 
sured at 77i^) to M§ predicted by QCD compared 
to the direct measurement. The evolution is car- 
ried out at 4 loops, while the flavour matching is 
accomplished at 3 loops at 2 m c and 2 mb thresh- 
olds. 



Figure 17. The ALEPH strange spectral function 
with various contributions indicated. Higher- 
mass contributions from final states with larger 
multiplicities lead to very small branching frac- 
tions and are estimated. 



6. Strange Spectral Function and Strange 
Quark Mass 

The spectral function for strange final states 
has been determined by ALEPH [5]. Fig. 17 
shows that it is dominated by the vector K* (890) 
and higher mass (mostly axial-vector) resonances, 
leading into a poorly defined continuum region, 
however in agreement with the quark model. 

The total rate for strange final states, using 
the complete ALEPH analyses supplemented by 
results from other experiments [35] is determined 
tobeB(T -> iy T hadrons s= -i) = (29.3±1.0) 1(T 3 , 
leading to 

i? T ,s = 0.163 ±0.006. (16) 

Spectral moments are again useful tools to un- 
ravel the different components of the inclusive 
rate. Since we are mostly interested in the spe- 
cific contributions from the us strange final state, 
it is useful to form the difference 

a kl — ^ pfci ^ jjkl 

- FT? \2 U T,S=0 - JTJ ]2 U r,S=-l { L <) 



where the flavour-independent perturbative part 
and gluon condensate cancel. Fig. 18 shows the 
interesting behaviour of A°° expressed differen- 
tially as a function of s. 

The leading QCD contribution to is a term 
proportional to the square of the strange quark 
mass at the r energy scale. Special attention must 
be devoted to the perturbative expansion of this 
term which has the early behaviour of an asymp- 
totic series. We quote here the recent result from 
the analysis of Ref. [36] , yielding 

m s {m 2 T ) = (120±ll cxp ±8 yila ± 19 th ) MeV(18) 

where the dominant uncertainty is from the- 
ory, mostly because of the poor convergence be- 
haviour. Stability checks have been performed 
and their effect has been conservatively taken into 
account in the estimate of the theoretical uncer- 
tainty. 
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Figure 18. Differential rate for A°°, differ- 
ence between properly normalized nonstrange 
and strange spectral functions (see text for de- 
tails). The contribution from massless pertur- 
bative QCD vanishes. To guide the eye, the 
solid line interpolates between bins of constant 
0.1 GeV 2 width. 



7. Conclusions 

The decays r — > v T + hadrons constitute a 
clean and powerful way to study hadronic physics 
up to y/s ~ 1.8 GeV. Beautiful resonance analy- 
ses have already been done, providing new insight 
into hadron dynamics. Probably the major sur- 
prise has been the fact that inclusive hadron pro- 
duction is well described by perturbative QCD 
with very small nonperturbative components at 
the r mass. Despite the fact that this low-energy 
region is dominated by resonance physics, meth- 
ods based on global quark-hadron duality work 
indeed very well. 

The new ALEPH preliminary results using the 
full LEP1 sample have been presented. Sat- 
isfactory agreement with CLEO is observed in 
the 7T7r° and 3irir° decay modes. The r spec- 
tral functions have now reached a precision level 
where detailed investigations are possible, partic- 
ularly in the most interesting tttt channel. The 
breaking of SU(2) symmetry can be directly de- 
termined through the comparison between e + e~ 
and t spectral functions. The difference between 
the two spectral functions does not agree with 



the computed effect of isospin symmetry break- 
ing and requires further experimentation. The 
determination of the widths of the charged and 
the neutral p is quite sensitive to the fit condi- 
tions, even though the data cover a very large 
mass range. 

The measurement of the vector and axial- 
vector spectral functions has provided the way 
for quantitative QCD analyses. These spectral 
functions are very well described in a global way 
by 0(al) perturbative QCD with small nonper- 
turbative components. Precise determinations of 
a s agree for both spectral functions and they also 
agree with all the other determinations from the 
Z width, the rate of Z to jets and deep inelastic 
lepton scattering. Indeed from r decays 

a s (M|) T = 0.1202 ± 0.0027 (19) 

in excellent agreement with the average from all 
other determinations [37] 

a s (M|)„ on _ r = 0.1187 ± 0.0020 (20) 

The strange spectral function yields a very 
competitive value for the strange quark mass 
which can be evolved to the usual comparison 
scale 

m s (2 GeV) = (116± 2 °) MeV (21) 

The use of the vector r spectral function and 
the QCD-based approach as tested in r decays 
improve the calculations of hadronic vacuum po- 
larization considerably. Significant results have 
been obtained for the running of a to the Z 
mass and the muon anomalous magnetic moment. 
Both of these quantities must be known with high 
precision as they have the potential to give access 
to new physics. 
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